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We propose experimental schemes to create and probe minimum forms of different topologically 
ordered states in a plaquette of an optical lattice: Resonating Valence Bond, Laughlin and string- 
net condensed states. We show how to create anyonic excitations on top of these liquids and 
detect their fractional statistics. In addition, we propose a way to design a plaquette ring-exchange 
interaction, the building block Hamiltonian of a lattice topological theory. Our preparation and 
detection schemes combine different techniques already demonstrated in experiments with atoms in 
optical superlattices. 
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Strong correlations between particles can lead to un- 
conventional states of matter that break the traditional 
paradigms of condensed matter physics P] . Among these 
exotic phases topological liquids [2] are at the frontier 
of current theoretical and experimental research. They 
are disordered states that do not break any symmetries 
when cooled to zero temperature. Surprisingly, they ex- 
hibit some kind of exotic order, dubbed topological order 
[3, [sj , which can not be understood in terms of a local 
order parameter. This global hidden pattern is revealed 
in the peculiar behavior both of the ground state, with a 
degeneracy that depends on the topology of the system 
and of the elementary excitations, which are anyons with 
fractional statistics [4]. 

The interest in topological liquids started in connec- 
tion with two landmark phenomena in condensed matter 
physics: fractional quantum Hall effect [5.] and high tem- 
perature superconductivity 0. In fractional quantum 
Hall systems electrons organize themselves in topological 
liquids, like the Laughlin state [7], following a global pat- 
tern that can not be locally destroyed. High temperature 
superconductivity was proposed by Anderson [8] to occur 
when doping a topological spin liquid: a Resonating Va- 
lence Bond (RVB) state in which the system fluctuates 
among many singlet bond configurations [8, 9*]. Recently, 
the study of topological states of matter has received spe- 
cial attention in the context of topological quantum com- 
putation [3, 10], which seeks to exploit them to encode 
and manipulate information in a manner which is resis- 
tant to errors. Moreover, understanding topological or- 
der may help us to understand the origin of elementary 
particles. According to Wen's theory [11], fundamen- 
tal particles, like photons and electrons, may be indeed 
collective excitations that emerge from a topologically 
ordered vacuum, a string- net condensate [ill]. 

Except for the fractional quantum Hall effect, there 
is no experimental evidence as to the existence of topo- 
logically ordered phases. It remains a huge challenge to 
develop theoretical techniques to look for topological liq- 
uids in realistic models and find them in the laboratory. 
In this direction, artificial design of topological states in 



the versatile and highly controllable atomic systems in 
optical lattices [12.1 ap pears to be a very promising pos- 
sibility 0, Q [la [li, ti3, H [13 . 

In this Letter, we show how to use ultracold atoms in 
optical lattices to create and detect different instances of 
topological order in the minimum non-trivial lattice sys- 
tem: four spins in a plaquette. Using a superlattice struc- 
ture lO, 111, [221, 0, 0] it is possible to devise an array of 
disconnected plaquettes, which can be controlled and de- 
tected in parallel. When the hopping amplitude between 
plaquette sites is very small, atoms are site localized and 
the physics is governed by the remaining spins. By com- 
bining different techniques we show how to prepare these 
spins in minimum versions of topical topological liquids: 
a RVB state, a Laughlin state, and a string- net condensed 
state. Making use of the ex per imental ability to control 
superexchange interactions jl^, [lH between neighboring 
spins a RVB state can be created [l^. This state can 
be transformed into a Lauhghlin state of two particles, 
which surprisingly appears in the absence of any rota- 
tion or effective magnetic field. This correspondence of 
the RVB state and the Laughlin state exactly demon- 
strates for the case of a plaquette the equivalence of 
these two states proposed by Laughlin [2^. To stabilize 
a string-net condensed state we develop a way to isolate 
a ring-exchange interaction [19] involving the four spins 
in the plaquette. The superlattice structure offers a set 
of precise tools to characterize these states. By merging 
plaquettes into double wells or single wells, in combina- 
tion with usual interference experiments and spin selec- 
tive measurements, the properties of these states can be 
revealed. Even though their minimum size, these states 
allow us to observe instances of topological order. By 
locally addressing each spin in a plaquette, we show how 
to create anyons on top of these liquids and detect their 
fractional statistics. 

As an additional application, a two-particle paired 
state with (i-wave symmetry, the fundamental compo- 
nent of a high temperature superconductor [6], could be 
created by doping a RVB state [14]. A novel technique 
allows us to reveal the characteristic d-wave symmetry of 
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FIG. 1: Schematics of optical lattice setup. Using an optical 
superlattice configuration along two orthogonal lattice direc- 
tions, an array of decoupled plaquettes can be created (a). 
The sites within a plaquette are denumbered as in (a). By 
controUing the two optical superlattices independently, differ- 
ent potential bias, and A^, can be introduced along the 
X— and y— direction, leading to different site energy offsets 
/Ji as well as different vibrational level splittings at the lattice 
sites. 



this state. 

The mini-topological liquids we consider here consti- 
tute fundamental building blocks of larger topologically 
ordered states 
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28|. Furthermore, plaquette ring ex- 
change interactions are the basic ingredients of lattice 
gauge theories [29], theoretical models describing topo- 
logical matter. By connecting plaquettes in the appropri- 
ate manner a variety of strongly correlated many body 
states could be achieved [27, 28]. 

The system. We consider a system of atoms in two in- 
ternal states (J =t, i, which for the case of e.g. ^^Rb 
atoms could correspond to the hyperfine states \F = 
l,mi? = +1) and \F = I^tuf = —1). The atoms are 
loaded into a two dimensional superlattice, which is pro- 
duced by superimposing a long and a short period lattice 
[22! both in the x— and direction in such a way that an 
array of disconnected plaquettes is created (see Fig. 1). 
In the following, we will restrict the discussion to the case 
of bosonic atoms, though similar results can be obtained 
in a straightforward manner for the case of fermions. 

The dynamics of atoms in a single plaquette is gov- 
erned by the Hubbard Hamiltonian 



H 



(id), 



where a^cr and riia- are, respectively, the bosonic annihila- 
tor and the particle number operator at site i and for spin 
a. By controlling the superlattice structure, the tunnel- 
ing amplitudes in the x- and y- direction, tx = ti2 = ts4 
and ty = ^23 = ^14, can be tuned independently. Fur- 
thermore, the dependence of the offset energies /i^cr on 
position and spin state can be designed using additional 
magnetic offsets or gradient fields. In the following we 
will make full use of the experimental ability to control 
these parameters, as already demonstrated in [22|,|23|] for 
a single double well. 



RVB in a plaquette. RVB states, in which particles 
are paired into short-range singlets, are one of the most 
relevant examples of topological spin liquids Q • We con- 
sider here minimum forms of RVB states consisting of 
four particles in a plaquette: 
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|0). 



(1) 



Here, slj = {a^^a^j^ — creates a singlet state on 

sites i and j and |0) is the vacuum state. The states 
([1]) are disordered states with zero local magnetization, 
(S^) = 0, for all sites i. They are both total singlets, 
with S-\^±) = Y.iS~\^±) = 0. But they behave dif- 
ferently under rotation of the plaquette by 90°: |^+) is 
even (has 5- wave symmetry), whereas |^_) is odd (has 
d-wave symmetry). As larger RVB states [sl, [30|, states 
([1]) exhibit topological order. 

In the following we develop a scheme to prepare and 
detect the state The state \^-) can be designed in 

a similar manner. We start with a situation in which we 
have four particles per plaquette and tunneling is only 
allowed along the y- direction. The system can be then 
prepared in a valence bond state |VB^) = SI4523IO), 
with singlets in the vertical bonds (see Fig. [2]). By adia- 
batically turning on tunneling along the x- direction we 
will connect the state iVB^y) to the state To make 

this connection possible, we consider a situation in which 
the tunneling amplitudes tx and ty are very small in com- 
parison to the on-site interaction energy U. Under these 
conditions, the particles are site localized and the physics 
is described by the superexchange interactions between 
the remaining spins: 



Hs = Jx (A,2 + ^3,4) + Jy (A,3 + ^,4) 



(2) 



Here, Vi 



is the projector onto a singlet state 



on sites i and j, Jx{y) = ^^e dots denote 

higher order terms in Jx{y)' M tx is suddenly increased 
from to ty the system will resonate with frequency 




FIG. 2: Preparation of a RVB state. Two spin-triplet atom 
pairs are split into double wells along the ^/-direction in a 
double well-potential (a). A magnetic field gradient along 
the ^/-direction is kept for a finite amount of time to turn 
the spin triplet bonds along the vertical direction into spin- 
singlet states (b-c). Finally, the tunnel coupling along the 
x-direction is increased adiabatically to connect the bonds, 
such that a RVB state is formed (d). The symbols for the 
bonds are explained in Fig. [3l 
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FIG. 3: (a) Eigenstates and eigenenergies of the spin Hamil- 
tonian ([2|) vs. tunnel coupling ratio along the horizontal vs. 
vertical bond, (b) Legend explaining the different symbols 
used to characterize spin singlet and spin triplet states of two 
particles. 

2Jy/h between the two valence bond states iVB^y) and 
|VBaj) = s\ 3IO). These oscillations could be detected 
by monitoring singlet pairs in the x- and y-bonds using 
a bandmapping technique [lH, 0] after merging the dou- 
ble wells in either the x- or ^/-direction into a single well. 
For the case of bosons, an initial spin singlet state with 
an antisymmetric spin wavefunction will lead to half of 
the population in the first excited state, whereas for a 
spin triplet state only the lowest vibrational state will 
be occupied after merging, thus allowing one to distin- 
guish between the two spin states [25] . In order to create 
the state |^+), tx has to be turned on adiabatically. As 
shown in Fig. [3l the states |VB^) and \^-\-) are adiabat- 
ically connected. The only state to which a transition is 
not forbidden by symmetry constraints is |^-), the other 
total singlet. The energy gap to this state is always on 
the order of ~ 2Jy^ giving a time scale of tens of ms, 
which can be easily fulfilled in experiments. 

The two RVB states presented in eq. ([T]), can be char- 
acterized and differentiated from each other using two 
alternative methods. 

a) Merging technique. By merging wells along 
the y- or diagonal direction we can monitor sin- 
glets in that direction via the bandmapping technique 
outlined in [22]. For example, by merging along 
the diagonals we could easily discriminate |<l>+) = 
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FIG. 4: Momentum distribution of the polarized two-particle 
states resulting from the RVB states |$+) (a) and |$-) (b) 
(see text). Such momentum distributions could be observed 
after releasing the atoms from the optical lattice potentials 
and a subsequent time-of- flight period. 

+ ^r,3^2^4 ~ 4,34,4) 1^)5 which is made out of 

triplets in the diagonal bonds, with t^j = 

from |$_) = s|3S2 4|0), made out of singlets. Alter- 
natively, we can distinguish these two states by merg- 
ing along the vertical direction. To do this we first 
undo the adiabatic path we followed before (see Fig. [3]) 
by decreasing tx from t to 0. In this way |^+) 
will be connected to the state 5|2'543|0), with singlets 
in the vertical bonds, whereas \^-) will connect to 

75 (^1^2^4,3 + ^M^ts ~ 4,24,3) with triplets in the 
vertical bonds. 

b) Conversion into a polarized two-particle state. Any 
four spin state |^) with well defined = J^- = can 
be written as a state of two up particles in a background 
of spin down particles 

|$)= ^ V(^i,X2)5+5+|UU>, (3) 

where is the spin raising operator on site x = 1, . . . , 4, 
and liiii) = a|^a2|a3^a4jO). If we remove the back- 
ground of spin down particles, that is, if we apply the 
operator Xli^j state ([3]), we are left with 

a system of two polarized hard-core bosons with wave 
function ?/^(xi,X2). In practice, spin down particles can 
be effectively removed by projecting the spin down part 
of state (j3j) onto the state Xl^^^ <^J|<^]jO|), where |0|) 
is the vacuum of down particles. The properties of the 
resulting two particle state directly reflect those of the 
spin parent state ([3]). By observing the momentum dis- 
tribution of the two particles via a common time of flight 
experiment, we can read back the spin-spin correlations 
of the parent spin state. For example, observing a hole 
at the center of the time of flight picture would be an un- 
ambiguous signature of a total spin singlet parent state. 
To see this, note that a spin state |$) is a total sin- 
glet if and only if (<l>|6'+6'" |<l>) = 0, which for the corre- 
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FIG. 5: Legend of correlated quantum states on a plaquette 
that are created using the methods described in the text, (a) 
s-wave RVB state, |$+), (b) d-wave RVB state, |$-), (c) 
LaughUn state, (d) paired state with d-wave symmetry. 



spending two-particle state translates into a zero occupa- 
tion of the momentum state with k = (0,0). Moreover, 
for the case of a total singlet, the spin structure fac- 
tor, a E.-.-e^^-^"""'^(45^.-), (with k = f(4,^^), 
X = {ix^iy)a^ ^x-,^y = 0,1), is directly equivalent to 
the momentum distribution of the corresponding two- 
particle state, since for a rotationally invariant state we 
have {S^S^') oc Yie{SiS^,). Additionally, the symme- 
try of the spin state (s-wave, p-wave, or d-wave) is in- 
herited by the two particle wave function and is di- 
rectly reflected in the center of mass angular momen- 
tum (Lcm = 0, ±1,2) of the pair. To see this, note 
that the center of mass angular momentum operator, 
^cm oc X]^=i i^ equivalent to the translation op- 

erator by one site along the ring, which rotates the pla- 
quette by 90°. 

This conversion technique provides us with a way to 
fully characterize the states ([1]). The corresponding time 
of flight pictures (see Fig. [4j) show both a dip at the 
center, probing their being total singlets, but are quite 
different in structure, a consequence of the different sym- 
metry of the 5- wave and d-wave RVB states. The two 
particle states obtained from the RVB states ([1]) turn 
out to be of special importance by themselves. We show 
below that the state \^^) can be converted into a Laugh- 
lin state, whereas \^-) can be transformed into a paired 
state with d-wave symmetry. 

Laughlin State. The Laughlin state \?\ is one of the best 
known examples of a topologically ordered state. Here, 
we describe a way to prepare a Laughlin state of two 
particles starting with the state |<l>+) in ([T]). In contrast 
to previous schemes our method does not require either 
rotation [sij of the optical lattice or of the individual 
wells [32^ , or the presence of effective magnetic flelds [is] . 
The key point is to realize that the state |<l>+) can be 
written in the form ([3|) with 



spin down particles, a Laughlin state of the remaining 
spin up particles is created (|4j). This state is an eigen- 
state of the total angular momentum operator, L = 
E^maJ^a^, with eigenvalue L = 4. Here, the opera- 
tor aj^ = I Y^t=i e^^^^^^aj, creates a particle in a state 
of angular momentum m. It has also a well deflned cen- 
ter of mass angular momentum, Lcm = 0. This vanishing 
Lcm is a consequence of the s-wave symmetry of the wave 
function, inherited from the state |^+). The state (|4]) is 
indeed a Laughlin quasihole state |7] . It contains a quasi- 
hole at the center of the plaquette, whose characteristic 
density proflle, with a dip at the center, could be observed 
in a time of flight interference experiment (see Fig. [3|). 

The equivalence of a long-range RVB state of 2N spins 
and a Laughlin state of N hard-core bosons has been pro- 
posed by Laughlin for a triangular two dimensional lat- 
tice Hamiltonian [2^. Indeed, it is known that this con- 
nection is exact for a lattice of spins sitting on a ring and 
interacting with a long-range interaction, the so called 
Haldane-Shastry model [ssj]: i^HS = j Jij^i.ji with 
J-~^ = sin \_i^{xi — Xj)\ . A surprising fact in our case is 
that the Laughlin state (|4]) appears in the absence of any 
frustration or long-range interaction (there is no inter- 
action between spins along the diagonals in Hamiltonian 
([2])). The key point to understand why the Laughlin 
state appears nevertheless, is to realize that for the spe- 
cial case of a plaquette the Hamiltonian ([2]) can be writ- 
ten in the form i^s = — ^HS, with S being the total 
spin operator. Since i^HS is rotationally invariant, we 
have that [i^HS, S'^] = and Hs and Hus share the same 
eigenstates with well deflned 5*^ . This has important con- 
sequences for the elementary excitations of Hamiltonian 
([2j). As for i^HS, they are ^-quasiholes spanned by wave 
functions of the form t/j^^ = {zi — r]){z2 — describing 
half of a boson missing at position r]. These quasiholes 
are ^-anyons according to the generalized deflnition of 
fractional statistics introduced by Haldane Q . In terms 
of spins, they are ^-spin excitations, the so called spinous 

Spinous and fractional statistics. Low-energy exci- 
tations of RVB states are created by breaking one of 
the spin-singlet bonds. For the case of a plaquette, 
they are = 1 excitations of the form 1^771,772) = 



V^(xi,X2) = ZiZ2{zi - Z2Y 



(4) 



where zi = e*2^% Xi = 1,...,4. By removing the 



|0), containing a pair of spinous localized 
at sites r]i and r]2- If the position 771 of one of the 
spinous is flxed, there are three possible states of this 
form, corresponding to three different positions of the 
other spinon (Fig. [6]). It is interesting to observe that 
these states are not linearly independent. They generate 
a subspace of dimension two which is orthogonal to the 
state Ei=r72,a.i,a;2 ^^"1 TTTT)- This non-orthogonality of 
states describing spinous at different positions is a char- 
acteristic feature of quasiparticles with fractional statis- 
tics, as deflned by Haldane In contrast to the case of 
spinous, states describing bosons or fermions at different 



(a) 





r 


\3 




(c) 






\72 







(c) 



FIG. 6: Spinon Excitations. Linearly dependent states con- 
taining two spinon excitations (marked by the two up arrows) 
on a plaquette. 



positions would be linearly independent. 

The fractional character of spinons becomes more 
transparent by mapping the spin system into a hard-core 
boson problem, in the same way that we did above. We 
consider the triangle obtained by excluding the site r]i 
in which one of the spinons is fixed. By removing down 
particles in this triangle, the spinon state is mapped onto 
a two-particle state of the form: 

i^m ^ ^1^2 (^-^ - f]2^ (^-^ - f]2^ {zi - ^2)^ (5) 

where Zi = e*^^% and Xi = 1,2,3 enumerates the sites 
of the triangle in consecutive order. The state ([5]) de- 
scribes a fractional ^-Laughlin quasiparticle [7] located 
at position 7^2- Since the addition of a complete boson is 
equivalent to a spin flip and creates an excitation with 
5*^ = 1, the quasiparticle, which constitutes half a boson 
corresponds to a spinon, with S = 1/2. 

Let us design an experiment to create spinons and 
probe their fractional statistics by detecting the non- 
orthogonality of the states 1^771,772)- ^ state of this form 
can be prepared experimentally by starting with the 
valence bond state sl^^^^sl.^^^JO), created as explained 
above. If the spin of the particle at site r]i is flipped this 

state is transformed into (^l^^^al^^^ + ^mi^mi) 4i,x2lO)' 
which by measuring the spin at rji in the z basis can be 
finally transformed into 1*^*771,772)- 

To probe the linear dependence of these three states, 
we project the down part of each of them onto the state 
^•ajjO), a projection that will yield zero in all cases 
only for linearly dependent states. 

Paired state with d-wave symmetry. Cuprate supercon- 
ductors are known to exhibit pairing with d-wave sym- 
metry [6]. A single pair with this exotic symmetry is 
described by the state: 



Ix) 
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'3,4 



4,4) |o). 



(6) 



Let us design a scheme to create and detect this state. 
Starting with the four-particle state |^-), we first re- 
move particles with spin down to obtain the state 
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^m) I^)' ^ triplet pair with the de- 



sired d-wdive symmetry. It is curious to see that the wave 
function describing this state, ex ^1^2(^1 + ^2)^(^1 — ^2)^, 
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FIG. 7: Eigenstates of the gauge lattice Hamiltonian (|7| (see 
text). Ground state (a). Flux quasiparticle excitation (b). 
Charge quasiparticle excitation, in which two bonds are ex- 
cited (c). 



corresponds to an excited Laughlin state, with a quasi- 
particle excitation in the center of the plaquette (the fac- 
tor z\Z2) and an excitation of the center of mass of two 
units of angular momentum (the factor (^1+^2)^)- This 
Lcm = 2 directly reflects the (i-wave symmetry of the 
state. In order to achieve the state (j6j) the triplet pair 
has to be transformed into a singlet. This can be don e by 
using the experimental techniques demonstrated in [36[. 
In order to reveal the (i-wave character of the state (|6]) we 
propose a novel technique which exploits the connection 
between the symmetry of the state and the center of mass 
angular momentum of the pair. By inverting the process 
above we transform the state ([6]) into a spin polarized 
pair. We then merge the four sites of the plaquette into 
a single well and convert the pair into a molecule using a 
photoassociation technique [37(. Since angular momen- 
tum of the center of mass is conserved in the merging 
process, the molecule will carry two units of angular mo- 
mentum, which will directly reflect the d-wave symmetry 
of the state ([6]). 

The paired state with d-wave symmetry (|6|) can be 
converted into a pair with a non- vanishing d-wave order 
parameter through e.g. admixture of a vacuum state. 
This could be done, by adiabtically increasing the tun- 
nelling to an unoccupied layer of empty plaquettes, below 
or above the occupied plaquette layer. 

Ring- exchange interactions. Lattice gauge theories [29| 
play an essential role in describing topological matter [31] . 
The minimum lattice gauge Hamiltonian describes a sys- 
tem of four spins in a plaquette and has the form: 



TT T CX QX QX QX i T \ ^ nZ Q Z 

Hg = -Ju ^1 ^2 *^3 + 2^ 



(7) 



It consists of four terms that commute with each other. 
The first one is a ring-exchange or flux interaction involv- 
ing the four spins. For < it favors symmetric states 
with respect to spin flipping of the whole plaquette. The 
other ones are charge interactions between neighboring 
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FIG. 9: By applying a magnetic field gradient along a diagonal 
direction of the plaquette, super exchange interactions can be 
suppressed. 



FIG. 8: Anyon Braiding. A flux- type and a charge- type quasi- 
particle are created by applying the operators S\ and 5*1 , 
respectively, to the ground state (a). A single charge- type 
quasiparticle is then moved around the flux quasiparticle by 
subsequent application of the local operators (b), S% (c), 
and SI (d). 



spins, which for J+ > favor states with anti-parallel 
neighboring spins. As for the case of an infinite lattice 
[lO[ the elementary excitations of this Hamiltonian are 
anyons. Though this is a well known result, for the sake 
of clarity of our discussion bellow, let us first briefly ex- 
plain it for the case of a single plaquette. The ground 
state of Hamiltonian ([7j) is a GHZ state of the form 



P) = ^(iTm) + imT)). 



(8) 



a maximally entangled state of four particles. It is in- 
deed the minimum version of a string- net condensate [llj , 
the ground state of ([7j) when extended to an infinite lat- 
tice. We can create two types of excitations on top of 
the state ([8]). They are flux-like or charge-like quasi- 
particles, (see Fig. [7|) depending on which term of the 
Hamiltonian ([7j) is excited. For example, a flux-like ex- 
citation (fluxon), which we denote by has the form 
!□) = -1= (Itltl) - litit)). It can be obtained by ap- 
plying, for example, the operator S\ to the state !□). 
Charge-like excitations, in which two neighboring spins 
become parallel, are always created in pairs. For exam- 
ple, the state !•□) = (|TiTT) + liTii)) contains two 
charge-like quasiparticles, one at the 1-4 bond and the 
other at the 3-4 bond. This state is obtained by apply- 
ing the operator S% to the state !□)• A charge-like quasi- 
particle can be moved around a flux-like one (see Fig. [8]) 
by applying the ring operator 8^828^8% onto the state 
!□). Since 5f5f5f5f|H) = 5f 5f 5f ^f^f p) = 
the wave function picks up a minus sign during the pro- 
cess. Therefore charges and fiuxons are relative ^-anyons 
in this model. 

In our optical plaquette a Hamiltonian like ([7j) seems, 
in principle, is difficult to implement. The reason behind 
is that four-spin interactions result from fourth order pro- 
cesses (higher order terms denoted by dots in equation 
([2])) in which four tunneling events occur. These pro- 



cesses are usually highly suppressed (~ /U^) compared 
to second order processes (~ /U)^ leadin g to dominant 
next neighbor superexchange interactions [l^, [HI- Here, 
we present a scheme to suppress second order processes 
in a plaquette, obtaining a dominating four-body inter- 
action. This will allow us to implement Hamiltonian ([7j) 
within a certain subspace of the spin Hilbert space. 
We consider a situation in which we have applied a mag- 
netic field gradient A along one of the diagonals (e.g., 
1-3) of the plaquette (see Fig. If A > At^/U, spin 
exchange interactions between neighboring sites are sup- 
pressed. The only remaining processes are either the ones 
in which the four spins in the plaquette are flipped (see 
Fig.) or those in which spins along the diagonal 2-4 are 
exchanged, giving rise to the Hamiltonian: 



Hr = -Jn (5+52-53+^4- +i7.c.)+J+^ 8^8] 

<'ij> 

+ Jx {8^8^ + H.c.) - A ^ 5,5f , 



(9) 



where Jn ^ 24^^//7^ J+ ^ 4tV/7, Jx ~ 16t'^//7^ and 

Bi = 0, ^2 = ^4 = 1, ^3 = 2. 

Within the subspace generated by the states |TiTi) and 
mil), Hamiltonian ([9]) is equivalent to Hamiltonian ([7j). 
The string- net condensed state !□) and the flux excita- 





FIG. 10: Effective ring exchange interaction in the presence 
of a magnetic field gradient. Fourth order processes connect- 
ing the states |||||) and |||||) are resonant (a-b). Second 
order processes connecting the states |||||) and |||||) are 
off resonance (c). 



tion !□) are therefore eigenstates of our system. Let us 
show how to prepare this states. We proceed as fohows. 
The plaquette is initiahy prepared in the state |TiTi). 
This can be done by starting with the state |TTTT), and 
then spin flipping atoms on the diagonal sites by address- 
ing them with the scheme presented in Fig. [Til In the 
presence of the magnetic field gradient, the system will 
evolve under Hamiltonian ([7j), oscillating between the 
states I Titi) and | itit) with a frequency uo = Jn/fi. 
For typical experimental parameters this frequency is of 
the order of 20 Hz, and can be resolved experimentally 
[22I . [2^. After an evolution time T = ir/Auj the system 
will be prepared in the state ^ (iTiti) + ^liTiT))^ a max- 
imally entangled state that can be easily transformed into 
either !□) or This can be done by applying the local 
phase operator Rq = e^^'^^^^^^\ with = 7r/4(-7r/4), 
which is performed by addressing sites 1 and 3, and let- 
ting the system evolve in the presence of a magnetic field 
B in the z direction for a time T = Oh/B. We can use 
Hamiltonian ([9]) together with local addressability of the 
plaquette sites to artificially create and detect the any- 
onic quasiparticles of Hamiltonian ([7j). Our proposal has 
the same spirit of the one recently proposed in [38[ , where 
anyonic states are artificially encoded using four photons. 
Even though anyonic states are as in [38] not eigenstates 
of our system the preparation and detection scheme we 
present here can be used in cases in which the Hamilto- 
nian ([7j) may be achieved using other methods [19]. Our 
scheme follows the idea proposed in ^] for anyon de- 
tection in small rotating atomic gases, a) Initialization. 
We prepare the system in the state (!□) — a su- 

perposition of a non-excited and a fiux-excited plaquette. 
Such superposition state results indeed from time evolu- 
tion of the state |TiTi) under Hamiltonian ([9]) after a time 
T = tt/Auj. b) Statistical phase accumulation. We then 
excite a pair of charge-like excitations and move one of 
them around the plaquette. This operation is performed 
by the operator Sf^f (see Fig. [8]), which we apply 
by subsequently addressing and acting on each site of the 
plaquette. Because of the relative —statistical phase of 
anyons, the state |H) will pick up a minus sign, and the 
system will end up in the state ^ (!□) c) Detec- 

tion. We finally let the system evolve under Hamiltonian 
([7j) for a time T = tt/Auj, obtaining the final state |TiTi). 
If the excitations happened to be bosons or fermions with 
trivial statistics, the final state would have been |TiTi). 
These two states can be easily discriminated by, for ex- 
ample, measuring Si. 

In conclusion, we have presented a collection of 
schemes to create and detect instances of topological or- 
der in a minimum system: a plaquette filled with two 
or four particles in an optical lattice potential. Many 
of these could be directly implemented in current exper- 
iments using the presently available manipulation and 
detection techniques. Furthermore, the plaquette Hamil- 
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FIG. 11: Single site addressing within a plaquette. By adjust- 
ing the superlattice potentials in the x- and ^-direction, the 
potential wells on the edges of a plaquette can obtain different 
vibrational splittings (a), e.g. huji 7^ h(jtJ2 7^ h(jtJ3. This can be 
exploited to target the spin on a single site and manipulate 
it without affecting the neighboring spins in the plaquette. 
In order to achieve this, Raman transitions (b) resonant to 
an intermediate excited vibrational state on a plaquette edge 
can be used. For sufficiently spectrally narrowband Raman 
pulses, the transitions will only be driven on a chosen single 
plaquette site. 



tonians we have considered could be used as unit opera- 
tions to, together with an increased optical resolution to 
resolve individual plaquettes, engineer large scale topo- 
logical liquids. 
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